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We propose a scheme for continuous detection of itinerant microwave photons in circuit quantum 
electrodynamics. In the proposed device, a superconducting qubit is coupled dispersively to two 
resonators: one is used to form an impedance-matched A system that deterministically captures 
incoming photons, and the other is used for continuous monitoring of the event. The present 
scheme enables efficient photon detection: for realistic system parameters, the detection efficiency 
reaches 0.9 with a bandwidth of about ten megahertz. 

PACS numbers: 42.50.Pq 03.67.Lx, 85.25.Cp 


Microwave quantum optics using superconducting 
qubits and transmission lines, which is realized in circuit- 
quantum-electrodynamics setups, is one of the hottest 
research area in modern quantum physics [ij. Exploit¬ 
ing the large dipole moment of superconducting qubits, 
circuit QED enables various quantum-optical phenomena 
that have not been reached by quantum optics in the vis¬ 
ible domain. In particular, we can readily construct opti¬ 
cal setups with excellent one-dimensionality , which 
are suitable to construct a scalable quantum circuit. 
However, the lack of an efficient detector for itinerant 
microwave photons has been a long-standing problem, 
and several approaches have been proposed to date. One 
approach is to capture a propagating photon determinis¬ 
tically into a resonator mode and detect it afterward. In 
recent experiments, the possibility of such capturing has 
been demonstrated with an excellent fidelity @ . How¬ 
ever, this approach requires precise temporal control of 
the system parameters that depends on the exact pulse 
shapes of the signal photons. Another approach is to use 
the Kerr effect mediated by superconducting qubits 0- 
[^, which may enable non-destructive photon detection. 
However, it has been revealed that a high distinguisha- 
bility of the signal photon number can be achieved only 
by cascading several identical qubits with negligible pho¬ 
ton loss in between, which is a challenging technical task 
presently [1,0. 

An alternative approach is to use the deterministic 
switchin g o f a A system induced by individual pho¬ 
tons [lOl - ll^ . which has been experimentally realized in 
one-dimensional systems Cllil- Note that this occurs 
as a result of single-photon dynamics: the destructive 
interference between the input and the elastically scat¬ 
tered photons enables the deterministic operation. Re¬ 
cently, detection of propagating microwave photons has 
been demonstrated using a A system realized in a tilted 
washboard potential of a current-biased Josephson junc¬ 
tion [iM3- A problem with this scheme could be the 
substantial dissipation upon detection and the resultant 
long dead time before resetting. More recently, we real- 
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FIG. 1: Schematic of the single-photon detector. A qubit is 
coupled dispersively to two resonators. Resonators A, B, and 
the qubit are respectively coupled to waveguides A, B, and C. 
We input a signal photon through waveguide A, a probe field 
through waveguide B, and a drive field through waveguide C. 


ized a A system formed by the dressed states of a qubit- 
resonator system and discussed its performance as a pho¬ 
ton detector [13, in, ill. This detector attains a high 
detection efficiency within the detection bandwidth, re¬ 
gardless of the signal pulse profile and with negligible 
dark counts. However, this detector should be operated 
in the time-gated mode, since the drive field to generate 
the A-type transition must be turned off during the qubit 
readout to obtain a high fidelity. 

In this study, we present a practical scheme for con¬ 
tinuous detection of itinerant microwave photons. We 
couple two resonators to a qubit: one is used for form¬ 
ing a A system Bin and the other is used for con¬ 
tinuous qubit monitoring (T3 - [^ . The proposed device 
enables continuous operation of the photon detector, pre¬ 
serving the advantages of our previous scheme [l^ , such 
as a high detection efficiency, insensitivity to the signal 
pulse shape, and short dead times after detection. More¬ 
over, the efficient detection is possible without cascading 
qubits [1,0- 

We consider a device in which a superconducting qubit 
is coupled to two resonators A and B (Fig. [T|). Setting 
h = 1, this system is described by TLsys = oJaCL^o. -\- 
ojbb^b -\- ujqa^a -\- ga(a^cr -f a^a) -|- gb{b^(J -\- a^b), where 
a, b, and a respectively denote the annihilation opera- 





















2 


tors for resonators A, B, and the qubit. Wq, <Db, and uiq 
are their bare frequencies, and ga and are the qubit- 
resonator couplings. We set ga, gb)/‘2'!T = 

(10,12, 5, 0.5, 0.4) GHz for concreteness. Since this sys¬ 
tem is in the dispersive regime, TLsys is rewritten as 

T-l-sys = (uJatt^a + uJbb'^b)aa^ + -|- (wa - ‘2Xa)a'^a 

+{ujb - 2xf,)5^5]cricr, (1) 

where Xr = 9rIi^r — bJq) is called the dispersive shift (r = 
0 , 6 ), and the renormalized frequencies of the resonators 
and the qubit are uJa = ijJa + Xa, + Xb, and uiq = 

^q-Xa-Xb- Their values are (xa,Xb)/27r = (50,23) MHz 
and {uja,bJb,bOq)/2 t: = (10.050,12.023,4.927) GHz. 

The qubit and the resonators are respectively cou¬ 
pled to waveguides, through which we apply three kinds 
of microwaves (Fig. [1}. Through waveguide C, we ap¬ 
ply a continuous drive held to the qubit to generate an 
“impedance-matched” A system by the dressed states of 
the qubit and resonator A. Through waveguide A, we 
input a signal photon to be detected, which determinis¬ 
tically induces a Raman transition and excites the qubit. 
Through waveguide B, we apply a continuous probe held 
for the dispersive readout of the qubit state. 

We denote the radiative decay rates of resonators A, 
B, and the qubit by Ka, Kb, and 7 ', respectively. Ka de¬ 
termines the bandwidth of the photon detector, which 
should be smaller than or comparable to the level sepa¬ 
ration of the dressed states [|3) and |4) of Fig. [2](b)]. Kb 
determines the phase shift of the probe held upon rehec- 
tion and Kb ~ 2xb is favorable for qubit readout [ 2 ^ . We 
set {Ka, Kb) 1271 = (20,46) MHz. Additionally, the qubit 
undergoes non-radiative decay and its overall decay rate 
7 often dominates 7 '. The photon detection efficiency is 
sensitive to 7 . We assume a reasonably long-lived qubit 
of 7/27r = 0.01 MHz (Ti Ri 16 gs) 

We engineer the dressed states of qubit-resonators sys¬ 
tem through the qubit drive. Theoretically, the qubit 
drive is described by 'Hdr{t) = + E‘^{i)o'], 

where Ed{t) = Ede~'“^‘^*^ is a monochromatic drive held. 
In the frame rotating at utd, we obtain a static Hamilto¬ 
nian, 

'Hsys+dr = {bOaO'a + Wb 6 '^ 6 )CTCT^ -|- [{uJq - UJd) 

+{LOa - 2xa)a^a + {uJb - 2xb)b'^b]a^a + Hd(tT'l' -|- a),{2) 

where fid = \fYEd is the Rabi frequency of the qubit 
drive. Hereafter, fid represents the drive power. 

First, we consider the case with fid = 0. The eigen¬ 
states of Hsys+dr are the Fock states \q,na,nb), where 
q{= g, e) denotes the qubit state and Ua and nb{= 
0,1, • • •) denote the resonator photon numbers. To hnd 
the optimal drive conditions, we restrict ourselves to the 
zero- and one-photon states [Fig.[2Ka)]. In this study, we 
use resonator A to form a A system and resonator B 
as a readout resonator that preserves the qubit state 
upon transitions. For this purpose, we should realize 
the level structure of Fig. ura), where a;|g_o,o) < W|e,o.o) < 
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FIG. 2: Dressed-state engineering of the qubit-resonators sys¬ 
tem. (a) Level structure of the bare states {fid = 0) in the 
rotating frame, (b) Level structure of the dressed states at 
the operation point {fid = fTJ^^). (c) Dependences of the 

decay rates on the drive power. The drive frequency is set 
at u}dl2n = 4.832 GHz. An impedance-matched A system 
is formed at ffj^^f2Tv — 10.75 MHz. {k^) is normalized 
by Ka {Kb), (d) Amplitude of the reflection coefficient |rs| of 
a continuous signal field applied through waveguide A, as a 
function of the signal frequency lJs and the drive power fid- 
The upper (lower) curve represents 0741 (oisi). 

<^|e,i,o) < W|g_i_o) (nesting regime for resonator A) and 
W|g,o,o) < W|e.o.o) ^|g, 0 .i> < W|e,o,i> (un-nesting regime 
for resonator B) [Q. This is done by setting the drive 
frequency within the range of ujq — 2xa < < ^q — 2xh- 

Next, we consider the case with fid > 0. The drive 
field mixes the bare states to form tffie dressed states. We 
denote the dressed states by |1), |2), • • • and label them 
from the lowest in energy [Fig.j^Kb)]. The states |1) and 

12) are made of the zero-photon states. From Eq. ([2|), 

they are given by 

|T) = cos6»i 2|5,0,0) - sin6»i2|e,0,0), ( 3 ) 

|2) = sin6li2|5,0,0)-bcos6»i2|e,0,0), (4) 

012 = arctan[ 20 d/(w|e,Q,o) - W|g,o,o))]/ 2 . ( 5 ) 

13) and |4) (|5) and | 6 )) are made of the one-photon 

states of resonator A (B), which are obtained by re¬ 
placing | 5 , 0 , 0) and |e,0, 0) in Eqs. ©-([SI) with |e, 1,0) 
and | 5 , 1 , 0 ) (Ig, 0,1) and |e, 0,1)). The radiative de¬ 
cay rate from \j) {j = 3,4) to |f) {i = 1,2) is given 
by kT = Ka\{j\a^\i)\^■ We confirm that = K 42 , 
K 32 = K 41 , and + K 32 = K 44 + k %2 = Ka- Namely, 
|3) and |4) decay in two directions, satisfying the sum 
rule of decay rates. Similarly, = Kq 2 , K 52 = 
and K 54 -I- K 52 = + Kq 2 = Figure ©c) plots kT 

and kF as functions of the drive power. In the drive-off 
limit {fid —t 0 ), K 32 = t Ka, 'k\i = k |2 Kb, and 
others vanish. This represents the simple decay of the 
resonator modes preserving the qubit state [Fig. ©a)]. 
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FIG. 3: Microwave response to a signal photon. The signal 
photon has a carrier frequency of liJbI2h = 10.05 GHz and a 
Gaussian pulse shape with the length of I = 100 ns. (a) Time 
evolution of the qubit excitation probability Pe{t). The probe 
power (uh) is indicated. The signal photon profile |/s(t)P is 
also shown in units of (8 In 2/7r/^)^^^. (b) Time-averaged qubit 
excitation probability p^{t) for (nt) = 0.05. The integration 
time At is indicated. 


As we increase the drive power, the decay rates for res¬ 
onator A are inverted, whereas those for resonator B re¬ 
main almost unchanged. This is because of our choice 
of the drive frequency ujd- At in Fig. [UJc), the 

four decay rates concerning resonator A become iden¬ 
tical. Then, a resonant signal photon deterministically 
induces a Raman transition of |1) —>■ \j) —>■ |2) {j = 3,4). 
Regarding resonator B, we should make K 52 (= as 
small as possible to suppress the |1) —>■ |fc) —>■ | 2 ) and 
|2) —>• \k) —)► |1) transitions (fe = 5,6). For this pur¬ 
pose, uJd close to ujq — 2xa is advantageous. We set 
Wdl^TT = 4.832 GHz [5 MHz above {ujg — 2xa)/27r] here¬ 
after, which results in H™^/27r = 10.75 MHz. Then, 
cos^ 0 i 2 = 0.99, cos ^034 ^0.61, and cos^^se = 0.96. 
Namely, the dressed states |1), |2), |5), and | 6 ) are almost 
identical to the bare states |(;, 0 , 0 ), |e, 0 , 0 ), \g,0, 1 ), and 
|e, 0,1), respectively. k| 2 (= ^ei) about 0.9% of Kb- 

In Fig. UDd), we plot the reflection coefficient jr^l of 
a weak continuous field of frequency ujs through waveg¬ 
uide A. We observe that impedance matching (|rs| ~ 0) 
takes place at fid — and W 31 ^ Ws < W 41 , where 

Wji denotes the transition frequency between \j) and |i). 
This indicates that each signal photon induces the Ra¬ 
man transition in the A system and is absorbed deter¬ 
ministically. We choose these drive power and signal fre¬ 
quency as the operating points of the photon detector. 

Next, we study the response of the detector to a single¬ 
photon signal. The signal photon is assumed to be a 
Gaussian pulse with length I and frequency ujg, namely, 
f,{t) = (81n2/7rP)i/4 x .^^hich is nor¬ 

malized as J dt\fs{t)\‘^ = 1. Setting u)sI‘2tt = 10.05 GHz 
and I = 100 ns, we plot the time evolution of the qubit 
excitation probability Pe{t), which represents the pop¬ 
ulation of 12), by the red solid line in Fig. EJa). Pe{t) 
increases within the pulse duration and approaches to 
unity, which agrees well with dt'\fs{t')\‘^. A high ef¬ 
ficiency is attained regardless of the pulse shape as long 
as the linewidth of the photon is narrower than that of 


the A system. After the pulse duration, Pe{t) decreases 
gradually by natural decay of the qubit with rate 7 . 

Now we consider the effect of the probe field, Ep{t) = 
. by which the qubit is continuously moni¬ 
tored |19l - [^ . From Eq. ([T|), the resonant frequency of 
resonator B depends on the qubit state. This is reflected 
in the phase shift of the probe field upon reflection. The 
phase shift 0g for the qubit ground state is given by 
0g = 2 arctan[Kb/ 2 (a;f, — Wp)], and 0e for the qubit ex¬ 
cited state is obtained by replacing^Wb with uJt — 2xb- 
Hereafter, in order to suppress the |1) —>■ |5) transition, 
we set the probe frequency at Wp = Wf, — 2xb rather than 
the usually chosen condition Wp = ujb — Xb- This results 
in e*®*' = (3 -I- 4i)/5 and e*®® = —1. We measure one 
quadrature of the reflected field for discriminating the 
qubit state. We infer the qubit state through the time- 
averaged probe field with an integration time At. The 
signal-to-noise ratio (SNR) is given, assuming that the 
noise is purely of quantum origin, by 


SNR = 


Kb{nb)At 



( 6 ) 


where (ub) = 4|Epp//tb represents the probe power in 
terms of the mean photon number in resonator B M- 
The readout fidelity E is given by 


E = erf(SNR/y2), (7) 

where erf denotes the error function In practice, the 
noise could be enhanced by technical reasons such as the 
noise from the amplifiers. Here, we assume the noiseless 
phase-sensitive amplification preserving the SNR (2^. 

In Fig. EKa), we plot the time evolution of Pe{t) in 
the presence of the probe field with various power (ub). 
We observe that the near-deterministic qubit excitation 
[red solid line in Fig. EKa)] is degraded by increasing the 
probe power. This is attributed mainly to the enhanced 
qubit decay through the | 2 ) —>• jfc) —>• | 1 ) transition (fc = 
5,6). However, for the probe power considered here, the 
backaction of the qubit readout is not severe and the 
qubit excitation is maintained for several microseconds. 
Hereafter, we fix the probe power at (ub) = 0.05. Then, 
E = 0.99 (0.999) [SNR=2.58 (3.29)] is attained by taking 
At =575 (939) ns. The long qubit lifetime enables us to 
take such long integration times. 

The single-photon detection efficiency 77 , which is the 
probability to find the qubit excitation, is evaluated as 
follows. Since we infer the qubit state through the time- 
averaged amplitude, we introduce the time-averaged ex¬ 
citation probability, 



and find the maximum probability p^(tm) in time. Con¬ 
sidering that the probability to correctly infer the qubit 
state is E) 12, p is given by 


,1-|-J^ . ,, 1 — 

V=Peitm)^r- + [^-Peitm)\ -- 


2 


2 


( 9 ) 
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FIG. 4: Single-photon detection efficiency. The probe power is fixed at (n^) = 0.05. (a) Dependence of the detection efficiency 
on the pulse length I for various integration times At. Values of the corresponding readout fidelity T are also indicated. 
The input photon is tuned at ldsI^-k = 10.05 GHz. (b) Detection efficiency as a function of Dd and tUs for I — 100 ns and 
At = 575 ns. The dashed lines indicate uiii, 0 J 31 , and (c) Cross section of (b) at I2tt = 10.75 MHz (red solid). 

The results for different drive conditions are also shown: uidl^n = 4.841 GHz and I2tv = 17.27 MHz (green dashed) and 
ajd/Svr = 4.850 GHz and j2n = 21.00 MHz (blue dotted). 


77 ~ 1/2 for a low SNR, implying that the qubit state is 
completely indistinguishable. In contrast, 77 ~ Pe{tm) for 
a high SNR. 

In Fig. [3Kb), we plot Pe(t) for various At. Expectedly, 
Peit) becomes flatter as we increase At, which implies the 
loss of detection signal. However, owing to the long qubit 
lifetime, the decrease of Pg(t) due to the time-averaging 
is at most several percent. In Fig. SKa), we plot the ef¬ 
ficiency as a function of the pulse length I of the signal 
photon. If the qubit lifetime is infinite, the efficiency in¬ 
creases monotonically with /. In practice, the efficiency 
is maximized at a finite pulse length due to the qubit 
decay during the pulse duration. We observe that a high 
efficiency is maintained for a wide range of /, which is an 
advantage of the A-based scheme. The loss of efficiency 
is due to the infidelity of the qubit measurement for short 
At, whereas it is due to the time-averaging for long At. 
For At = 575 ns {J- = 0.99), the maximum efficiency of 
0.91 is obtained at 1 ~ 90 ns. In Fig. SKb), we plot the 
efficiency as a function of fid and uig. Comparing this 
with Fig. |2^d), we confirm that the impedance-matching 
leads to a high detection efficiency. The cross section 
of Fig. HKb) at fid = is shown in Fig. UKc) by the 

red solid line, which shows the detection band of this de¬ 
tector. The detection efficiency exceeds 0.9 (0.8) for a 
bandwidth of 9 MHz (20 MHz). 

Four final comments are in order, (i) The detector is 
insensitive to signal photons when the qubit is excited. 
This causes a dead time of the detector, which amounts 
to several microseconds at {nt) = 0.05 [Fig. jSKa)]. How¬ 
ever, by applying a reset pulse upon detection of the qubit 
excitation |18j . we may shorten the dead time to several 
hundreds of nanoseconds, (ii) The detection band center 
is tunable by changing the drive condition. We show the 
detection band for different drive conditions in Fig. UKc). 
The detection band has two peaks located at W 31 and 
W 41 in general. However, as we increase uJd and accord¬ 
ingly , K 52 and kqi are increased. This enhances the 
probe backaction and degrades the detection efficiency, 
(iii) In the continuous measurement, one may worry that 
the quantum Zeno effect prohibits efficient photon detec¬ 


tion, since the apparent measurement time interval seems 
infinitely small. However, even in the continuous mea¬ 
surement, the effective measurement time interval Atm 
remains finite, which is determined by the dephasing rate 
induced by the measuring apparatus [^[^. Here, the 
probe field functions as the apparatus and Atm is deter¬ 
mined by SNR^l, namely. Atm ~ ^/Kb{nb) = 175 ns. 
This is obviously long enough to avoid the Zeno effect 
[see Fig. EKa)]. (iv) The probe photons may cause the 
dark counts by inducing the | 1 ) — \k) —>• | 2 ) transition 
(fc = 5,6). We can numerically check that this proba¬ 
bility is about 0.2% per one probe photon. Therefore, 
the dark count rate is estimated to be (142 ps)~^ for 
(ub) = 0.05. A lower dark count rate is accomplished by 
reducing the probe power. 

In summary, we proposed a practical scheme for con¬ 
tinuous detection of itinerant microwave photons. The 
detector consists of a qubit and two resonators in the 
dispersive regime. We apply a drive field to the qubit to 
form an impedance-matched A system, a signal photon to 
one of the resonators to excite the qubit, and a probe field 
to the other to continuously monitor the qubit. For real¬ 
istic system parameters, the detector has a maximum de¬ 
tection efficiency exceeding 0.9 and a bandwidth of about 
ten megahertz. One can improve the performance of the 
detector further by increasing the qubit lifetime and/or 
the dispersive shifts. 

This work was partly supported by MEXT KAKENHI 
(Grant Nos. 25400417 and 26220601), Project for Devel¬ 
oping Innovation Systems of MEXT, National Institute 
of Information and Communications Technology (NICT), 
and ImPACT Program of Council for Science, Technol¬ 
ogy and Innovation. 

Appendix A: Detection efficiency 

In the main part of this study, we defined the detection 
efficiency intuitively through the time-averaged qubit ex¬ 
citation probability. Here, we investigate the detection 
efficiency more rigorously on the basis of the quantum 
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jumps of the qubit. We observe that the deviation be¬ 
tween these two definitions is negligible for the parameter 
range discussed in this study. 

1. Time-independent case 


a. Method 1 

In the main part of this study, we evaluate the de¬ 
tection efficiency as follows. From the qubit excitation 
probability Pe{t), we define the time-averaged probabil¬ 
ity pe{t) by 


In the dispersive readout of the qubit state, we infer the 
qubit state through the time-averaged probe field. First, 
we preliminarily observe a case in which the qubit keeps 
staying in its ground/excited state. We denote the field 
operator for the probe port by c(t), which is normalized 
as [c(t),c’l’(t')] = 6{t — t'). We apply a classical field 
(coherent state) as the probe of the qubit state. The 
phase of the probe field is sensitive to the qubit state as 


m) 


EpC^^o for I 5 ) 
for |e) ’ 


(Al) 


where the natural phase factor e is neglected. We 
introduce the time-averaged field operator by 


c(t) 


I 


f 

Jt-At 




(A2) 


which is normalized as [c(t),c’^(f)] = 1. We infer the 
qubit state through one of its quadratures, x{t) = 
Im[e“®^®»“''®'=^/^c(t)], which maximizes the signal-to-noise 
ratio (SNR). The expectation value of this operator is 


{x{t)) 


—Epy/At sin I 5 ) 

Ep\fKi sin ^ ^ for |e) 


(A3) 


We set the threshold at {x) = 0 and judge the qubit 
state through the sign of {x). Since c{t) is normalized as 
[c(t),c^(t)] = 1, the quantum noise in each quadrature 
is 1/2 for a coherent state. The SNR and the readout 
fidelity are then given by 

SNR = 2Ep'/Ktsin(^-^^^'^ , (A4) 

T = erf(SNR/y2), (A5) 


which are Eqs. m and © of the main part. The proba¬ 
bility to correctly infer the qubit state is (I + E)I2. 


2. Time-dependent case 

Next we investigate a more realistic situation in which 
the qubit is excited at t ~ 0 and decays gradually af¬ 
terward [Fig. 3(a)]. The detection efficiency rj is defined 
as the probability to detect the qubit excitation. We 
compare two methods for evaluating this probability: In 
method I, which we adopted in the main part of this 
study, we intuitively evaluate the detection efficiency iji 
through the time-averaged qubit excitation probability. 
In method 2, we evaluate the detection efficiency 772 more 
rigorously based on the quantum jumps of the qubit ob¬ 
served in actual measurements. 


^ / cli'Pe(i'): (A6) 

Jt-At 

and find the moment tm that maximizes pe(t). We define 
the detection efficiency 771 as the probability to detect the 
qubit excitation at this moment. Considering that the 
probability to correctly infer the qubit state is (I-|-J^)/2, 
771 is given by 


T?1 = Peitm) 


l+E 

2 


[1 


Pe(tm)] 


1-E 

2 


(A7) 


b. Method 2 

The qubit excitation/de-excitation is observed as the 
quantum jumps in actual measurements. We consider a 
single event where the qubit is excited at t = 0 and is 
de-excited at t = r. Considering the rapid response time 
of the resonator (1 /k{, ~ 3.5 ns), we may regard that the 
probe field responds immediately to the quantum jumps 
of the qubit as 


m) 


Epe^^<’ (t < 0, T < t) 

EpE^‘ (0 < t < r) 


(AS) 


{x) is maximized at t = (r-l-At)/2. The maximum value 
depends on the duration r of the qubit excitation as 


{x) = 


- Epy/Kt sin {1-2^) (0<T<At) 

Ep^fKi sin ^ ^ {At < t) 


(A9) 

Accordingly, the probability q{T) to detect the qubit ex¬ 
citation is 


q{t) 


O' 

1 -erf 

SNR 7-, 

. C 2 

11 

1 -f erf 

snr'^ 
V2 )\ 


(0 < T < At) 
{At < t) 

(AlO) 


The shape of q{T) is shown in Fig.ISja). It is a monoton- 
ically increasing function of r and becomes constant for 
T > At. 

We denote the probability distribution function of the 
duration r of the qubit excitation by Q{t)^ which is nor¬ 
malized as dTQ{T) = I. Then, the overall probability 
772 to detect the qubit excitation is 


m 


dTQ{T)q{T). 


(All) 
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FIG. 5: (a) The probability g(r) to detect the qubit excitation 
as a function of the excitation duration time r (red solid), and 
its step-function approximation (blue dashed), (b) Compar¬ 
ison of rji (dotted) and ri 2 (solid). The qubit lifetime is 
assumed to be 3 fis (blue), 6 fj,s (red), and 16 /rs (green). 


We can show that rji and 772 are almost identical if 
the qubit decay within the integration time At is small, 
namely, TAt < 1. Regarding rji, pe{t) is maximized at 
tjn = At, and the maximum value Pe(im) is approximated 
well by Pe{At/2). Equation (IA7I) is then rewritten as 


771 =p,(At/2)l^ + [1 -p,(At/2)]l^. (A14) 


On the other hand, regarding 772 , we may replace g(r) 
by a step function [dashed line in Fig. [Sl[a)] as long as 
Q{t) is almost constant for 0 < r < At. Then, using 
lAt /2 d'TQ{T) = pe(At/2), 772 is rewritten as 


c. Comparison of rji and 772 

Here we compare 771 and 772 assuming a simple case in 
which the qubit is excited at t = 0 and decays exponen¬ 
tially with a rate of T. The excitation probability Pe{t) 
is given by 

Peit) = (A12) 

where 9{t) is the step function. The probability distri¬ 
bution Q{t) is connected to Pe{t) by dTQ{T) = Pe{t). 
Q{t) is therefore given by 

Q(r) = Te-^A (A13) 


772 =Pe(At/2)i^ + [1 -p,(At/2)]i^, (A15) 

which coincides with rji of Eq. (IA14I) . 

In Fig. EKb), 771 and 772 are plotted as functions of the 
integration time At, setting the probe power at (rib) = 
0.05. We confirm that the deviation between 771 and 772 
is small for At < T”^. From Fig. [3l[a), we estimate that 
the qubit lifetime is about 6 ps for (rib) = 0.05. Then, 
\Vi ~ ^ 2 ! ^ 0.01% for At < 1 ps. Thus, we can safely 
regard rji as the detection efficiency. 
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